Can oscillating scalar fields decay into particles with a large thermal mass? 
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We calculate the dissipation rate of a coherently oscillating scalar field in a thermal environment 
using nonequilibrium quantum field theory and apply it to the reheating stage after cosmic inflation. 
It is shown that the rate is nonvanishing even when particles coupled to the oscillating infiaton field 
have a larger thermal mass than it, and therefore the cosmic temperature can be much higher than 
inflaton's mass even in the absence of preheating. Its cosmological implications are also discussed. 
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In our contemporary understanding, the origin of the primeval fireball whose existence was assumed in the conven- 
tional hot big bang cosmology is the reheating processes after inflation — an accelerated cosmic expansion which has 
made the universe homogeneous and spatially flat with small density fluctuations that eventually grow to the observed 
large-scale structure 0, Q ■ The universe is reheated through the dissipation of coherent oscillation of the zero-mode of 
the infiaton, the scalar field whose potential energy drives inflation. While the initial stage of reheating could be rather 
complicated due to an explosive particle production induced by parametric resonance, which is dubbed as preheating 
1^], the final stage is dominated by perturbative decay. The latter process determines the reheat temperature, Tr, 
the temperature at the outset of the radiation domination |jj . 

Note, however, that in general Tfj is much lower than the highest temperature the universe has ever experienced 
after inflation even in the case only perturbative decay operates to reheat the universe @. This means that in the 
late stage of the reheating processes, the inflaton decays not in a vacuum but in a thermal medium. About this point 
an interesting claim has been made in |^ that if the would-be decay products of the oscillating inflaton acquire a 
thermal mass larger than the inflaton mass in the thermal background, it cannot decay into these particles, and that 
reheating is suspended for some time, based on the observation that the phase space would be closed for the mass of 
the decay product being larger than half the inflaton mass. The decay width of the inflaton (j) with mass into two 
massive particles with mass m reads 



where F^o is the decay rate in the case m = 0. So if we simply replace m with a thermal mass m{T) ~ gT and if it 
is larger than m^/i, the phase space is closed and inflaton decay is apparently forbidden. Here g is some coupling 
constant of the would-be decay product. Then thermal history after inflation would be drastically changed. That is, the 
highest temperature in this era cannot exceed ~ w^/g if preheating is inoperative, and also the reheat temperature is 
bounded from above by / g and is independent of the decay rate of the inflaton in case the conventional calculation 
gives a larger value. The former would change the abundance of supermassive particles and the latter affects the 
gravitino abundance |^ , because the gravitino-entropy ratio after inflation is proportional to the temperature at the 
onset of radiation domination. 

Furthermore, this situation is not specific to the reheating stage after inflation but may apply in any epoch when 
significant amount of entropy is produced out of the decay of oscillating scalar field with a relatively small mass. 
Indeed the above possibility was first pointed out by Linde 7] in the context of Affleck-Dine baryogenesis Q where 
the Affleck-Dine scalar field oscillates with a mass of order of 10^^"^ GeV in a medium with a much higher temperature. 
The final magnitude of baryon asymmetry changes if this suspension of decay is operative J7j . 

The above naive picture, however, may be too simplistic because a thermal mass is different from the intrinsic mass 
and because coherent field oscillation is different from a collection of particles. Thus it is very important to analyze 
this problem from a more fundamental point of view, for it has a profound implication not only to the cosmology of 
the early Universe but also to particle physics in that it strongly affects various species of the particles produced in the 
early universe as mentioned above. In this paper, extending our previous work Q, we analyze this problem in terms 
of a nonequilibrium quantum field theory at finite temperature [Tol |. Inclusion of thermal masses of the would-be 
decay products of the infiaton is achieved by adopting a resummed propagator when we calculate the effective action 
for (j). As a result of resummation the self energy of the decay product acquires not only real part, which appears as 



I. INTRODUCTION 
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a high-temperature correction to the mass, but also an imaginary part. The latter plays a crucial role in determining 
the dissipation rate of the inflaton. Consequently, we find that the inflaton can dissipate its energy even when its 
would-be decay products have a larger thermal mass than the inflaton itself. 



II. MODEL AND EQUATION OF MOTION 



For clarity we adopt a simpler model than 9J, that is, we adopt the following Lagrangian. 
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Here </> is an oscillating real scalar field and x is another real scalar field. can decay into a pair of x particles 
through the interaction A^(/)x^, if it is energetically allowed. The dimensionful coupling constant M. may be written 
as = hm^ in some supersymmetric inflation models where /i is a Yukawa coupling . In such models can also 
decay into two fermions, which is suppressed by Pauli-blocking at finite temperature. On the contrary, the decay into 
two bosons is enhanced due to the induced emission. This is the reason we consider the latter decay process. 

We analyze the behavior of the above system under the following assumption which mimic the cosmological situ- 
ations we are interested in, such as the late reheating phase after inflation. First we neglect cosmic expansion since 
we are interested in the phenomena which occur in a shorter time scale than the expansion time. Second we assume 
X is in a thermal state with a specific temperature T = and that it acquires a large thermal mass due to the self 
coupling. Note that x can easily be thermalized during the reheating stage since its thermalization rate, ~ g'^T, can 
naturally be much larger than the cosmic expansion rate. Finally the scalar field (/> is oscillating but we consider the 
situation the parametric resonance is already terminated with a field amplitude A^|0| < vr?^. 

Due to its coherent nature, scalar field oscillation behaves nearly classically, but its decay is of course a quantum 
process. So we calculate an effective action for </) and derive an equation of motion for its expectation value. For 
this purpose we should use the in-in or the closed time-path formalisms in which the time contour starting from the 
infinite past must run to the infinite future without fixing the final condition and come back to the infinite past again 
in calculating the gene rating functional 12]. This method has been applied to various cosmological problems by a 
number of authors jl.'j ll4L llq . The generating functional in the present model is given by 



Z[J, K\ = Tr 



T_ <^ exp 



i / dt (fx{J^(l>^ +K^x-] 



T+ < cxp 



i / dt \d^x[J+(j)+ + K+X+) 



JW[J,K] 
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where denotes a field component X on the plus-branch (— oo to +oo) and X^ that on the minus-branch (-|-oo 
to — oo). The symbol T+ represents the ordinary time ordering, and T_ the anti-time ordering. J± and K± are the 
external fields for and Xi respectively, p is the initial density matrix which is assigned according to the assumption 
above mentioned. 

In terms of the components along the plus and the minus branches, the effective action reads 



T[(j)+,(l)_] = W[J+,J-,K±^0]~ I dt (d^x[J+{x)(l)+{x) - J-{x)(j)-{x) 

J — OO J 

with (j)+{x) = 5W\J+,J-]/5J+{x) and (t>^{x) = -5W[J+, J-]/ 5J-{x). 



(4) 




FIG. 1: One-loop Feynman diagram incorporated in the effective action. Solid line denotes 0, and broken line x- 

Here we consider a one-loop correction depicted in Fig. ^ which includes the essential effect in our analysis for 
illustration. We also note that instead of 0+ and it is more convenient to use (j)c = -I- (/)_)/2 and 0a = 0+ — 4>- 
and set 0a — ^ in the end because 0+ and 0_ should be identified with each other eventually. Then the effective 
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action to this order is given by 

r[0c, 4>a] = - J d*x(j)A{x){0 + TO^ )0c(a;) 

d'^xd^x'Cix - x')9{t - t')(t)Aix)<j)cix') + ^ / d^xd'^x'D{x - x')<j)A{x)(t>Aix') , (5) 



Cix - x') = AM'^lm [G^ix ~ x')^] , D{x - x') = 2M''Re [G^{x - .t')'] , (6) 
where G^{x) is the Feynman propagator at finite temperature. Its Fourier modes read Il7| 

G^(p,t) = Jd^xG^{x,t)e-'-- = ^{[l + nB{u;p)]e-'^-'+nB{u;p)e'^-'}, tUp = -^p^ + (7) 



Up 

i 

— + ie 



G^ip) = dt G^ip, t)e'^°* ^ , _ ^ „ ^ + 27rn^K) - m^) , n^K) = (8) 



The resultant effective action is complex-valued as a manifestation of the dissipative nature of the system. We 
cannot obtain any sensible equation of motion by simply differentiating with respect to a field variable because we 
are dealing with a real scalar field and its equation of motion should be real-valued. The cure for this problem has 
been proposed by Morikawa [l^ . according which wc introduce an auxiliary random Gaussian field, ^{x), to rewrite 
the effective action as follows. 

exp(zr[0,>A]) = y"pe^'[^]exp{ireff[</>e,0A,C]}, (9) 

where 

T,s[(j)„(j)A,i]=Rer[(j)c,(l)A]+ f d*x^{x)c^A{x). (10) 



Here P[^] is a statistical distribution functional for ^(x) which is a Gaussian with its dispersion given by the imaginary 
part of the effective action, {^{x)(,{x') ) — D{x — x'). 

While the mathematical equivalence between the above decomposition and the original expression (jSJ can easily be 
confirmed by performing the path integral with respect to ^(a;), we keep it as it is to obtain a real- valued equation of 
motion through 



0(PA 

From ((Tn|l . it reads 



= 0. (11) 



{U + mD (j)c{x) + dt' d^x'C{x - x')(j)c{x') = £.{x) . (12) 

J —oo J 

Hereafter we omit the suffix c. The solution to the above equation of motion can be readily found through Fourier 
transform, 



0,(0 =y d'x<^(a;,^)e-"=■^ 4>u{uo) ^ j dtcp,{t)e'^\ (13) 
etc. We find p2|) is transformed as 

+ + ^ S^iu)] 4>^{u) - itoVkkH = e.(^), (14) 



where we have defined 



hico) ^ S.iu.) ^ f ^P^^C^c.'), (15) 

2Ll> I Ztt uj ~ u 
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with 



Ck[Lo)= / dtd^xC{x)e'' 



(16) 



which is pure imaginary. Since S'fc(tj) is divergent, we subtract the divergence at = to renormahze mass |l5j |. 
We can show that the renormahzed part S^rcni^) = Skii^j) — 3^(0) is of order of Ai'^ / (47r^). From now on we refer to 
as the renormahzed mass, and assume that Ai < racj^. Then one can neglect Su{i^) in (|14|l to yield the solution 



(j)k{U) cos Mk{t -ti) + 



sinMkit-U) 



e 2 



iffc(Mfc)(t-ii) 



+ ii i dt'e-i^^^''^^^'-'"> sinA4(i - i')e.(i'), Mfc ^ + fe^ 



(17) 



Here we have also assumed f ^.(Mj.) ^ M^, which is justified whenever perturbation theory applies. From the above 
solution we can read off that r/j(Mfc) gives the dissipation rate of the field oscillation. In particular, the dissipation 
rate of the zero-mode condensate is given by 



fo(™0) 



2m^ 

m0 



-,1/2 



1 + 2nB ( ^ 



rB(o) 



1 + 2nB 



2TTm\ 



2m^ 
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(18) 



(19) 



(20) 

which arises in turn because the spatial Fourier mode of the Feynman propagator at finite temperature can be 
expressed as (0). 

III. INCLUSION OF THE THERMAL MASS OF THE DECAY PRODUCTS 



for 777,0 > 2777^. Here the last approximate equality applies at high temperature T 3> m^, and we have used 

/d^p 1 
-—— — (1 + 277„) [5{m^ - 2u)p) - S{m^ + 2ujp)] , 77„ ee 71b(cj„), 

for fc = mode. Here the delta functions appear as a result of time-integral of the form 

COs(7770 =F 2u!p)dt — 27r (5(7770 =F 2Wp), 



So far, although we have taken into account the fact that the decay product x is thermally populated, we have not 
considered effects of the thermal environment on x itself which induces a thermal mass to x particles. If we simply 
replaced the intrinsic mass 777^ with the thermal mass m^{T) in H18(l and H19(l . we would find the dissipation rate 
would vanish for m-^(T) > m^/2 as was claimed in [^Q- Here we carefully examine how to incorporate the thermal 
mass of X in our analysis in order to see the validity of such a simplistic argument. To this end, we should not rely 
on the formula (|18|l but go back to the more fundamental equation (|19|) to reconsider how this equation was derived. 

First we note that x's mass in Up of (|19|l is that appears in the denominator of the Feynman propagator ^ and 
the delta function emerges due to the infinitesimally small imaginary part ie. The effect of the thermal environment 
on which gives rise to finite-temperature correction to its mass, can be incorporated to the calculation of the 
dissipation rate of (f) if we apply resummation and use a resummed propagator of x instead of the finite-tem per ature 
bare Feynman propagator to calculate the effective action. With the help of the Matsubara formalism fl^, the 
denominator of the propagator acquires a self energy whose real part yields a finite-temperature correction to the 
mass as desired. This resummation procedure, however, generates appreciable magnitude of the imaginary part to the 
self energy at the same time, so that the delta function seen in H19|) will no longer be present as we see below. Thus 
we expect that the use of a resummed propagator changes the result qualitatively, and that the simple observation 
that a large thermal mass would close the phase space of the decay rate of the inflaton lfTH|l would not apply. 

We write real and imaginary parts of x's self energy, E(p), as ^r{p) and E/(p), respectively. Then the spectral 
function reads 



1 



1 



(llJ + leY — — 777^ — E/f — i'Sj {uj — ie)^ — — 777?, — Ejj + zE/ 



(21) 
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where 



= p^ + ml + ^nip) + Tl^ito) ^p^+ ml + Efl(p), T^p{u;) = (22) 
The full dressed propagator is then given by 

{[1 + nei^)] e{t) + nB{^)e{-t)} p,{p,u:)e-'^* 

In 

- ^ { [1 + nsiJ^ ~ ^^^,)] e— ^1*1-^-1*1 + usi^J^ + ^^e^->\-^ } , (23) 

which should be compared with the bare propagator |(7J) T^. Then since the propagator has a complex phase now, 
the simple cosine integral in our previous calculation, which gave rise to the delta function in H19fl . is replaced by 

dt e-^^xP* cos(m, - 2^;)t = „ (24) 

namely, the Breit-Wigncr function, when we calculate IjlQI) in terms of the dressed propagator (|23() . 

Since (|24|) is finite even for ^ 2uj'p we find the dissipation rate of (j> is nonvanishing even when <C 2m^{T). 
To the first order in T^p, the zero mode dissipation rate in such a regime is given by 

^ ^ f d^p If 2r^p(l + 2n;) 2r^p(l + 2n; 



2m^ J (2^)3 Lo'p^ \ ATIp + (m^ - 2c^p2 4r2^ + (^^ + 2a;p2 



4r2p + m2 4r2p + (m^ - 2^;)2 4r2p + (m^ + 2^^ 



/ ^2 



(25) 



where rt^ = ns^LUp). All the terms in the second line vanish in the limit T^p — > 0, while those in the first line reduce 
to delta functions in the same limit and the previous result with the undressed propagator is recovered except for the 
replacement ujp — > uj'p. 

As mentioned above, the above result remain finite even if is smaller than 2ijj'p. In particular, if thermal mass 
of X is much larger than the inflaton mass, to' ^ m^^ we obtain 



ro(m,) = 4-, / dc.;^^-- (1 + 2n'p)T,p + / dco'p^ n'^in'^ + l)T,p. (26) 



In the present simple model, x thermalizes only through self interaction then T^p{ujp) is given by T^p{Lu'p) = 

Sg'^T^ / {128nu}'p) We find the second term yields the dominant contribution in (|26|l for -C m^{T) « gT/2. 

As a result we obtain 

" 2^24^- 

Here the former factor is identical to the dissipation rate to massless particles at high temperature T ^ m^, and 
the latter factor represents the suppression due to the large thermal mass of the decay product. One may wonder 
the suppression factor might be proportional to 5^ just as T^p^ujp). However, the result of integration of the second 
term in H26|l yields (T ) in the denominator, which partially cancels to g^. Although the numerical value of the 
suppression factor in (|27|l is specific to the present model, it is a generic feature that the suppression is proportional 
to some combination of coupling constants which is related to thermalization processes of the decay product, because 
it arises from T^p/uj'p. 



IV. DISCUSSION 



The resultant reheat temperature in the present model is given by 

^ / 90 Y^' M^Mg 3g2 _ n2/g. V^/y M \\ 

^""^[-^J ^^24^'^-°''^° M200J [WGcvJ [l¥Gev) ^^^^ 
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where 5* is the effective number of relativistic degree of freedom and Mq is the reduced Planck scale. Again it is 
suppressed by the same factor, 3(7^/(247r^), compared with the case the scalar field decays into massless particles 
in the high-temperature medium, see eq. (157) of 9]. The above result has been obtained under the assumption of 
7710 ^ qTr and M < m^. The former condition reads 

„,«l,„.l„»,(|_)-"(_^)""GeV^ ,29) 

The above reheat temperature (|28|l should also be compared with Tr « rn^jg which would apply in the case large 
thermal mass could forbid inflaton decay completely. 

Our result has important implications to the abundances of supermassive dark matter particles and gravitinos as 
well as baryon asymmetry. First, in case a large thermal mass forbids inflaton decay completely, the temperature 
after preheating cannot be higher than ^ m^/g. Then the abundance of supermassive particles with mass mx is 
exponentially suppressed as cx (mx / mcj,Y e~'^"^^ /^'i' for mx ^ and 3 = 10. Our result shows that such a 
suppression is absent and an appreciable amount of supermassive particles could be created after inflation. Second, 
the gravitino abundance is suppressed by a factor of Sg^/ (247r^) in the present model compared with the conventional 
reheating scenario, because its abundance is proportional to the reheat temperature. Hence we can relax constraint 
imposed by the gravitino decay to this extent. 

In summary, we have calculated the dissipation rate of an oscillating scalar field in a thermal bath such as the 
inflaton in the late reheating stage, and shown that it is nonvanishing even if the would-be decay products have a 
thermal mass larger than the mass of the oscillating field. This yields several important implications that has not 
been taken into account so far, as discussed above. 
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